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Abstract. Consider the problem to explicitly calculate the law of the first 
passage time T(a) of a general Levy process Z above a positive level a. In 
this paper it is shown that the law of T(a) can be approximated arbitrarily 
closely by the laws of T n (a), the corresponding first passages time for X n , 
where (X n ) n is a sequence of Levy processes whose positive jumps follow a 
phase-type distribution. Subsequently, explicit expressions are derived for the 
laws of T n (a) and the upward ladder process of X n . The derivation is based 
on an embedding of X n into a class of Markov additive processes and on the 
solution of the fundamental (matrix) Wiener-Hopf factorisation for this class. 
This Wiener-Hopf factorisation can be computed explicitly by solving iteratively 
a certain fixed point equation. It is shown that, typically, this iteration converges 
geometrically fast. 
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1 Introduction 

A Levy process is a stochastic process with cadlag paths and homogeneous 
independent increments. Levy processes have turned up in a number of areas 
of applied probability, for instance as models for the workload of queues, the 
content of a dam and the reserve level of an insurance company, see e.g. [3113111%] 
and references therein. More recently, Levy processes have also been proposed as 
models for the evolution of the logarithms of interest rates or prices of assets and 
as model for credit derivatives, see e.g. El ESI to name just three articles. 
In many of the mentioned areas one is interested in an (explicit) characterisation 
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of the distributions of the running maximum, the first passage time over a level 
and the ladder height and ladder time processes. 

The study of pathwise extrema of a stochastic process is called fluctuation 
theory. For Levy processes with only jumps of one sign, the fluctuation theory 
is simplified and results are more explicit than in the general case; See e.g. the 
comprehensive review of Bingham [§], Chapter VII in Bertoin [5] or the reviews 
from martingale ^3] or potential-theoretic perspective In the case of a 

general Levy process with positive and negative jumps few explicit results are 
available. However, restricting one-self to the class of Levy process with positive 
jumps of phase-type (see Section 12. II below for definition), Asmussen et al. |H] 
explicitly characterised the joint law of the first exit time from the negative 
half- line and the corresponding overshoot. Using this result, the law of the first 
passage time of a positive level for a general Levy process can be approximated 
arbitrarily closely, as shown in Section 12.21 Indeed, any Levy process can be 
approximated arbitrarily closely by a Levy process with phase-type jumps and 
this convergence is shown to carry over to the first exit times mentioned. 

In [S] martingale techniques were invoked to characterise the joint law of 
the aforementioned stopping time and overshoot as solution of a certain linear 
system. In this paper, we study the same class of Levy processes X but follow 
a different approach from that of p]. First, the Levy process is embedded into 
a class of Markov additive processes (or a phase-process perturbed by a spec- 
trally negative Levy process) and then the matrix Wiener-Hopf factorisation 
is characterised. This factorisation extends existing results in the literature: 
We mention in particular Asmussen [2] and Rogers [201, who covered the fac- 
torisation of random walk and phase processes perturbed by Brownian motion, 
respectively. 

The analysis of this Wiener-Hopf factorisation leads to an alternative char- 
acterisation of the law of the maximum of X and to an explicit description 
of the law of the upcrossing ladder process of X. For explicit calculation of 
this factorisation, an algorithm is derived and it is shown that, typically, this 
algorithm converges geometrically fast. 

The rest of the paper is organised as follows. In Section|21the notation is set 
and some theory regarding the Wiener-Hopf factorisation and ladder processes 
of Levy processes is revisited. Section [3] is then devoted to the matrix Wiener- 
Hopf factorisation of the process X and its embedding. In Sections 01 explicit 
identities for the law of the first passage time over a positive level and the law 
of the ladder process are derived. Finally, in Section an algorithm for the 
computation of the aforementioned laws and its convergence are studied. 

2 Preliminaries 

2.1 Phase-type distributions 

A distribution F on (0, oo) is phase-type if it is the distribution of the absorbtion 
time C in a finite state continuous time Markov process J = {Jt}t>o with one 
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state A absorbing and the remaining ones 1, . . . , m transient. That is, Fit) — 
P(C < t) where C = mf{s > : J s = A}. The parameters are m, the restriction 
T of the full intensity matrix to the m transient states and the initial probability 
(row) vector a = {a\ . . . a m ) where on = P( Jo = «)• For any i = 1, . . . , m, let t\ 
be the intensity of a transition i — > A and write t=(ti ... t m )' for the (column) 
vector of such intensities. Note that t = — Tl, where 1 denotes a column vector 
of ones. It follows that the cumulative distribution F is given by: 

1 - F(x) = ae Tx l, (1) 

the density is f(x) = ae Tx t and the Laplace transform is given by F[s] = 
f™e- sx F{dx) = a(sI-T)-H. Note that F[s] can be extended to the complex 
plane except at a finite number of poles (the eigenvalues of T) . Throughout we 
will assume that the representation of the form for the distribution function 
F is minimal, that is, there exists no number k < m, k- vector b and k x fc-matrix 
G such that 1 - F{x) = be Gx l. 

2.2 Phase— type Levy processes 

Consider a Levy process X of the form 

+*<->, (2) 

where = {X^\t > 0} is a Levy process without positive jumps and 

= {x[ + \t > 0} is an independent compound Poisson process with in- 
tensity \( + ) and jumps of phase-type with parameters (to^ , T^ + ' , a^). We 
exclude the case that X(~> is a negative deterministic drift. We define by 
k(s) = kx(s) = loglE[e sX ] the Levy exponent of X which is well defined at 
least for s with Si(s) = and which is in this case given by 

where n X {-) denotes the Levy exponent of X^' and -F( + '[,s] = aS + '{sI — 
j<(+)^—if i s Laplace transform of F^ . Note that k can be analytically ex- 
tended to the positive half-plane except a finite number of poles (the eigenvalues 
of — T^ +t> ) and we shall denote this extension also by k. 

Any Levy process L may be approximated arbitrarily closely in law by a 
sequence (X(n)) n of Levy processes of type (J2J. Indeed, it is well known that 
X(n) weakly converges to L (as a process, in the Skorokhod topology) if and 
only if (X(n)i) n converges to L 1 in distribution (see e.g. Jacod and Shiryaev 
[T3] . Cor. VII. 3. 6). Also, the set of phase-type distributions is dense in the 
set of probability distributions on (0, oo) (in the sense of weak convergence of 
probability distributions). Therefore the aforementioned approximation can be 
obtained by first restricting the Levy measure v of Lj to K\[— e, e] and then us- 
ing that the probability distributions c+l{2; >c }^(da:) and c~l{ x< ^ t }v(dx) (with 
cf = l/u(e, oo) and c~ — l/v(—oc,—e) and 1a the indicator of the set A) 
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can be approximated arbitrarily closely by phase-type distributions. The rele- 
vant methodology for fitting a phase-type distribution to a given set of data is 
developed in [7] for traditional maximum likelihood and in ^U] in a Bayesian 
setting. 

Write T + (a) — T + (a)(X) for the first passage time of X over a 

T+(a) = in£{t >0:X t >a} (3) 

and denote by + (a)(X) = X T +^ — a the corresponding overshoot of X. The 
next result shows that the weak convergence of the processes X(n) carries over 
to the stopping times T + (a)(X(n)) and the overshoots + (a)(X(n)) 

Proposition 1 Let (X(n)) n be a sequence of Levy processes such that X{n) — > 
X weakly in the Skorokhod topology. Then, as n — > oo, 

(T+(a)(X(n)),0+(a)(X(n)) — (T+(a)(X),0 + (a)(X)) 

where the convergence is in distribution. 

Before giving the proof we first review Wiener- Hopf factorisation of Levy pro- 
cesses. 

2.3 Wiener-Hopf factorisations 

In this subsection, we review some of the fluctuation theory of Levy processes. 
For more background we refer the reader to Bingham |Uj or Bertoin |S] . Denote 
by X( + > = {i : ^St(pi) > 0} the set of roots pi with positive real part of the 
Cramer-Lundberg equation 

K{p) = K X {p) = a. (4) 

Let x A = min{x, 0} and x V = max{x, 0} and write M t = sup s<t (W t V 0) 
and L t — in£ s <t{X t A 0) for the supremum and infimum of X up to time t, 
respectively. Similarly, write M.[ ~^ and if' for the the corresponding quantities 
of X^ . Denote by e(a) an independent exponential random variable with mean 
a' 1 . Set for s with K(s) > 

<Pa ( s ) = E[exp(s7 e(o) )], tpH-8) = E[exp(-sM e(a) )]. (5) 

The functions s i— > <p^(s) are analytic for s with ±3t(s) > 0, respectively. By 
bounded convergence it follows that 95^(00) = P(J e ( a ) = 0) and fa(~ 00 ) = 
P(M e ( Q ) = 0). For a > 0, the functions s 1— > (p%(s) satisfy the Wiener-Hopf 
factorisation 

a/(a-K(s)) =(p+(s)ip-(s) for all s with 3?(s) = (6) 

For a proof we refer to e.g. Bingham [SJ Thm. la]. Since |y„ (s)| = \ip~ (s)| < 1 
for s with 3?(s) = 0, there are no roots of J2J with zero real part when a > 0. If 
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X is a Levy process of the form |J5J, the phase-type representation is minimal 
and a > 0, then tp+ is explicitly given by (as shown in [5]) 



det(-£j-rW) a e x( +) (-^) 
det(-T (+) ) ' IligK+)(*-ft)' 



(7) 



where the first factor is to be taken equal to 1 if X has no negative jumps. 
The following Wiener-Hopf identity (e.g. Thm. le]) links the joint distri- 
bution of the first time of crossing the level a and the corresponding overshoot 
(T+(o), 0+{a)) to the Wiener-Hopf factor ip+: 



We have made all preparations for the proof of Proposition ^ 

Proof of Proposition Q] Writing M™ = sup n<s<t (X(n) s V 0) for the running 
supremum of X(n), the triangle inequality implies that |M" — M t \ is smaller 
than sup s<t \X(n) s V — X s V 0| which is smaller than sup s<t \X(n) s — X s \. 
Since, by assumption, X(n) converges to X in the Skorokhod topology and 
since a Levy process is continuous at each fixed times a.s., it can be directly 
verified (as in Prop. VI. 2. 4 in ^3|) that M" converges in distribution to M t 
for fixed t > 0. This implies that the moment-generating function (mgf) of 
MZ „) converges to the mgf of M e ^ q y By the Wiener-Hopf identity |S| and the ex- 
tended continuity theorem it follows then that the joint Laplace transforms 
of (T + (a)(X(n)),0 + (a)(X(n))) converges to the joint Laplace transform of 
(T + (a)(X),0 + (a)(X)), which shows the stated convergence in distribution. 
□ 

Closely related to the supremum process M are the ascending ladder time 
and the ladder height processes, which are Levy processes. One of the objectives 
of this note is to explicitly find the distribution of these two processes. To be able 
to describe the ladder process, we first need look at the local time of M — X 
at zero. We will distinguish between two different cases (see [SJ Ch. IV] for 
details): 

If X^~> is not the negative of a subordinator, is regular for M — X and the 
local time process can be taken to be continuous. The canonical choice in this 
case is to take L = M c , the continuous part of the supremum process M. 

If X^ is the negative of a subordinator, is irregular for M ~ X (i.e. 
Po(<7 = 0) = 0, where a is the first time M — X reaches 0) and the zero set 
of M — X forms a discrete set. The supremum is a jump process, where the 
jump sizes form an i.i.d. sequence and the jump times are precisely this zero 
set. In order to let it fit in the same framework as the previous case, an extra 
randomisation is needed, which can be thought of as the analog of the random 
time change to transform a random walk into a compound Poisson process. Let 
ei(c), 62(0), 63(c), ... be a sequence of i.i.d. exponential random variables with a 
certain parameter c > and denote by n(t) the number of zeros of M — X up to 
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yj-(-A) - 

9>g (-**). 



(8) 
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but excluding time t, i.e. n(t) — max{i : Ri < t} where Ri are the subsequent 
zeros of the process M — X. Then the local time L t at time t is given by 



n(t) 

L t = Y.en(c). (9) 

n=l 

The ladder time is the right-inverse L _1 of the local time L 

Lt 1 = inf{s > : L s > t} 

and is a stopping time with respect to the standard filtration generated by X. 
The ladder height process H is taken to be H = M L -i, the supremum at this 
stopping time. Write k + for the joint characteristic exponent of the ladder 
process (L~ l , H) 

exp{-K+(a, /?)} = IE [expj-alT^l) - 0H (1)}] , a,/3 > 0. 

Denote by Gt = sup{w < t : X u = M u } the last time before t that X was at 
its supremum and let as before e(q) be an independent exponential time with 
mean Then an extension of the basic Wiener-Hopf factorisation (see e.g. 
8, Ch VI] or 0) tells us that (G e{q) , M e(q) ) and {e(q) - G e{q) , M e{q) - X e{q) ) are 
independent and the joint distribution of (G e i q \, M e i q \) is expressed in terms of 
k + by 

E r e -aG e(5 , -bM <q) i = a , 6 > 0. (10) 

K+{q + a,b) 

Note that this expression and the stated independence imply that 

[ J " K +(q,0)q + a-K(bY [ ' 

Similarly, we can define kT to be the joint Levy exponent of (L~ l ,H), the 
ladder process of the dual X = —X of X and then we have the following relation 
between the different characteristic exponents k, k + , cp+ for s with K(s) = and 
a > (e.g. Bertoin 8, Ch. VI]): 

, , , K + (a,0) 

^ («) = + f v (12) 
«; + (a, — sj 

K (s) = -a-«; + (0,-s)«;"(0,s), (13) 

where C is some constant. The second identity in l|13|l is often referred to as 
the Wiener-Hopf factorisation of the Levy exponent k. 

2.4 Embedding 

To study properties of the Levy process X, we follow ideas of [S] and embed 
X into a Markov process (A,Y). Informally, ones get A from X by 'levelling 
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out' the positive jumps into piecewise linear parts of gradient +1; the process 
Y is set equal to zero if X moves like X^' and equal to the underlying phase 
process of an upward jump otherwise. 

More precisely and slightly more generally, let Y be an irreducible Markov 
process with finite state space E, with E — {0, 1, . . . , m}. Denote the generator 
of Y restricted to E by Q. Letting be the spectrally negative Levy process 
of J5J and v and a functions that map E to R, the additive functional A — 
{A t } t >o is defined as 

A t = A + [ a(Y s )dX^ + f v(Y s )ds. (14) 
Jo Jo 

The pure fluid model corresponds to setting a equal to zero whereas this model 
reduces to a spectrally one-sided Levy process if v = 0, a is constant and Y is 
recurrent. 

Next, wc define Y + to be the upcrossing ladder process of Y, that is, 

F+ = F( 7t ) where lt = inf{s > : A s > t}. (15) 

It is easily verified that this time-changed process Y is again a Markov process; 
we denote its generator by Q + . The next section is devoted to a characterisation 
of its form. 

Consider now the special case of above additive functional (|14H where m — 
roW, and restriction of the intensity matrix of Y to {0, 1, . . . , m} is given by 
(in block notation) Q , where for a > 

Qa=[ t (+) T ( + ) J- (16) 

Moreover, for i G E set v(i) = 1 — a(i) and a(i) = Soi, the Kronecker delta, and 
let To(t) = J Q I(Y S — 0)ds denote local time of Y at 0, that is, the amount of 
time before time t that Y has spent in state 0. Then we get back the original 
process X by time-changing A with the inverse local time 

T Q - 1 (u)=irf{t>Q:T (t)>u}. (17) 

Indeed, the independence of the increments of X implies that X is in law equal 
to A o Tq 1 . Exponential killing of the original Levy process X at rate a can be 
incorporated by replacing Q by Q a , since all states of Y other than originate 
from positive jumps of X so that the local time T of A at zero corresponds to 
the time scale of X. 



3 Matrix Wiener-Hopf factorisations 

Denoting by v and a functions that map E to R and [0, oo) respectively, we 
consider now the additive functional A = {At,t > 0} given by H14JI. In the 
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sequel, we restrict ourselves^ to a characterisation of the generator Q + of the 
upcrossing ladder process Y + of A given in (|14() . In the analysis, we shall 
distinguish between the cases that — X^ is a subordinator, (that is, X^' has 
non- increasing paths) or that X^ has non-monotone paths. We partition the 
state space E of Y into a part E~~ , where v(j) < 0, and a part where 
v(j) > 0. It follows that A decreases as long as Y is in one of the states in E~ 
and that A can decreases as well as increase when Y is E + and the state space 
of the upcrossing ladder process Y + defined in l|15fl is given by E + . The matrix 
of "up-crossing phase probabilities" is denoted by r), that is, 

r / (i,i)=F[y + =j\Y = i], ieE~, j£E+. (18) 

Finally, let K(a, — G) be the matrix whose rows are given by (liK X (-) (—a(i)G)) 
with Kjf(-) given by (by the Levy-Khintchinc formula) 

Kjrc-) (-G) = ^G 2 -cG+ [ (e- Gx -I- xGl {N<1} ) v{dx), (19) 

with v the Levy measure and s the Gaussian coefficient of '. 

Theorem 1 (i) If —X 1 -^ is a subordinator, the generator matrix Q + of the 
process Y + satisfies the matrix equation given, in block notation, by 

Q K(a, -G) +q(J)=V Q G, (20) 

(ii) If Q is transient or if Q is recurrent and sup t A t — oo a.s., the solution in 
G e Q(\E+\) of HJ) is unique. 

Now we turn to the case that A is given by (|14p. where — X^ - ' is not a 
subordinator. We restrict ourselves to the case that the function v is positive, 
v(i) > , for each i £ E for which a(i) = 0, to ensure that A can increase in 
each state in E. 

Theorem 2 (i) If —X^> is not a subordinator, the generator matrix Q + of 
the process Y + defined in f_?5l) solves the equation 

K(a,-G)+Q = VG. (21) 

(ii) If Q is transient or if Q is recurrent and suptAt — +oo a.s., Q + is the 
unique G £ <2(|-E|) that solves 

Proof of Theorem (i) Define the function / : E x M — * R by f(j, x) — 
Ej jX [/i(Y fc + )l{ 7fc<00 }] where h is any function on E and where Ej, x denotes 
the expectation under the measure P conditioned on {Yq = j,A = x}. Since 
Y + is a Markov process with generator Q + , the function / can be expressed as 

f(i,x) = l' i exp(Q + (k-x))h, i £ E, x < k. (22) 
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However, on the other hand, it is straightforward to check that /(Y" tA7fc , A tAlk ) 
is a martingale and we find by ito's lemma that / = (/(«, u),i E E) satisfies for 
i G E and u < k, 

T(a(i)f(i, u)) + v(i)f(i, u) + Y, Hi (fti, «) - /(*, «)) = 0, (23) 

j 

where T denotes the infinitesimal generator of the process X^~>: 

2 poo 

r/(^) = y/"(z) + c/'(z) + j (/(* + y) - /(*) - y/'( 2 ;)l { |,| <1} ) u(dy). 

Substituting equation ()22)l into equation 1)23)1 and using that h is arbitrary, we 
conclude that Q + solves equation J2U- (ii) Now we turn to the proof of the 
uniqueness of the solution of l|21l) . To that end, let G E Q(|-B|) be another 
solution of l|21|) and define the function / by replacing Q + by G in 1)22)1: by an 

application of Ito's lemma it follows then that f(Yt,At) is a local martingale 
that is bounded on [0,7^] and invoking the optional stopping theorem yields 
that 

f{j,x) =E i)X [/(y tA7fe!J 4 tA7 J] 

= E ijSB [/(?+ , A 7fe )l {7fc <oo}] + km E A B [J(F t , A t )l {7fe=oo} ] . (24) 

By definition of / and the absence of positive jumps of A, the first expectation 
in (|24|l is equal to f(j, x). Note that second term in l|24|l is zero if Q is transient 
or Q is recurrent and sup t A t — +oo. Indeed, in the latter case, jk is finite a.s., 
whereas in the former case P(Yf E E) converges to zero. Thus f = f and the 
statement (ii) follows. □ 
Proof of Theorem Q] Denote by the matrix Q + E Q(|i? + |) the generator of F + 
and let h be any function on E + . Then the Markov property of Y + yields that 
f:ExR^R defined by /(», a:) = E^ x [/(F+)l {7fc<oo} ] is given by 



l-exp(Q+(fc-a;))ft, 
l^exp(Q + (fc-a;))h, i E E' 



Following the line of reasoning of the proof of Theorem [21 and replacing every- 
where (|22|) by 125) 1 . we arrive at Theorem ^ □ 



3.1 Matrix factorisation of X 

In case (A, Y) is the embedding as in Section 12.41 of a Levy process X of the 
form (J2J) we can say more about the structure of the generator matrix Q + of the 
associated time changed Markov process Y + in 1)15)1 . We distinguish between 
the cases in which X^' has non-increasing or non-monotone paths. 

In the former case — X^' is a subordinator and the additive process A does 
not increase on the set {s : Y s — 0} and we have E~ = {0}. Note that in 
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this case the matrix T] in l|18|) is of size 1 X (\E\ — 1). To show explicitly the 
discounting at rate a we write 

Va(j) =E , [e^l {f + =j} ] ! a>0, (26) 

where IE^o denotes the expectation conditioned on {Aq = Yq = 0}. As before, 
for an additive process A of the form (|14|l to be the embedding of a Levy 
process X of the form J2J), we set a(i) — 8oi and v(i) = 1 — a(i) in (|14|l and let 
the intensity matrix Q a of the special form (|16|l . Inserting these quantities in 
equation 1|20|) with X = (<5oi, « G -E)diag and V = I — X leads to the following 
result: 

Proposition 2 // — is a subordinator, the generator matrix of the 

Markov chain Y + in }15\) corresponding to the embedding is given by 

Q+ = T ^+t^ Va , (27) 
where the vector r\ a satisfies the equation 

Va = #)qW ((A<+> + a)I k x( -> (-Q+)) _1 . (28) 

There is an alternative probabilistic derivation of the formula l|27|l for the inten- 
sity matrix which uses the Markovian structure. See [2] for a similar argument 
in a random walk setting. Let M(i,j) denote the ijth element of a matrix M. 
By the form of the time change, Y + is equal to Y when A is at a maximum. 
Thus, for i ^ j, the total intensity Q^ihj) of a jump of Y + from i to j, is 
given by the intensity T + (i,j) for the direct transitions from i — * j added to 
(t^ + ^r/J l )(i, j), the intensity of first passing from i to A and then being renewed 
with initial distribution r) a . The latter is equal to the distribution of Y + at the 
end of a negative excursion of A away from its supremum. Similarly, the total 
rate Q^(i,i) of a jump in state i is equal to the rate T + (z,z) of a direct jump 
added to the rate (t^t]' a )(i, i) of a jump after being killed and renewed. 

In the next section we shall also give a probabilistic argument to prove the 
second relation l|28|l . Now we turn to the case X'~) does not have decreasing 
paths. Let &(q) denote the largest real root of n X (-) (s) ~ q, which is positive 
for q > 0. 

Proposition 3 If —X^~^ is not a subordinator, the generator matrix of the 
Markov chain Y + in corresponding to the embedding is given by 

Q+ = M (+) +m (+) r7 a (29) 

where 

M(+) = + A<«) ^ ro(+) = _ M < + > 1= (*(°Y <+, >) (30) 
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and rj a satisfies 

\(+) 

Va = ^jTxw(°. a(+) ^x<-) Htf )• ( 31 ) 
where ip~ (— Q) = f °° e^P^ ^ € dx) /or generator matrices Q 

Proof The form of the matrices £ and V and equation l|21|) imply that all rows 
of from the second one on till the last one are given (in block notation) by 
(t(+) ) T (+) ). Rewriting J2H) we find for the first row 



A(+) 



("*<-)(-<#)- (« + A«)l) 



A« 

-0,a (32) 



A(- 



where Afi denotes the first row of a matrix M. From (|12|l with T^ + ' = 0, it 
<q) = snp s< e (q) 



follows that Mf I = sup s<e( >, Xs \ the supremum of X^ ) at an independent 



exponential time e(q), has an exponential distribution with mean <I>(g) , so 
that for q > and Sft(s) > 

9 _1 (^<-) (s) - g)vx(-) (*) = - s), 

where y^ ( _)(s) = f~ X (-) (s), the moment generating function of the infimum 
/^rj = inf s < e ( g ) xi ' of Af(~) at e(q). By the Cayley-Hamilton theorem this 
relation remains valid with s replaced by a non negative definite matrix (and 
thus in particular with s replaced by — Q„ ). Multiplying both sides of 13 2|) from 
the right with the matrix ^r(-) (— Qa) yields that 

\(+) 

*(a + A< + ))- 1 ($(a + AW)I + Q+)i - ^^(O, «K +A c + ),xC-) 

Thus we find that the first row of is given by 

\(+) 

(Q+)i - -$(a + A«) + *(a + AW)-^ (0, a)^ A(+)jAr( _) (-Q+) 

and the result - |J3TJ follows. □ 
As above there is also a probabilistic derivation of the form of . Write 
T\ for the last ascending ladder time of X^ before the first jump of X and Ti 
for the first ascending ladder time of X after T\ . and let the vector r] a denote 
the distribution of Y + at the end of the excursion away from the supremum 

r, a (fc)=E[e- Q ^- T ^l {yr2=fc} ]. 

Note that the supremum of X before the first positive jump of X (where X is 
killed at rate a if a > 0) has the same distribution as that of X^> killed at 
an rate a + A^ + '. Since Af(~) has an exponential with mean l(a + A(+)), the 
intensity Q+(0,i), i ^ 0, of a jump of Y + from — > i is the intensity <I>(a + A( + ') 
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of the jump — » A times the initial distribution 77+ (i) for Y + to be renewed 
in state i at the end of the excursion of A away from its supremum. As before 
we see that the rate Q+(0,0) of Y + jumping in state is the sum of the rate 
$(a + A (+) ) for Y + to jump directly and the rate $(a + A(+))?7 a (0) to exit first 
and be renewed to state 0. 

The other rows follow by the argument given in the subordinator case. In 
the next section we shall also give an alternative derivation of (|31l) . 

4 First passage and ladder processes 

If X is a Levy process of the form and X^ is not a subordinator, it follows 
from Proposition|21that the running supremum M e ^ = sup 0<t<e( - (; - l X t of X at 
an independent exponential time has a phase-type distribution given by 

P(M e(g) > k) = 1 exp { (M (+) + mWr/,)fc} 1, (33) 

where lo denotes the row-vector with a one in the position corresponding to 
E° = {0} and else. Using alternative proofs, this result was found before in 
Asmussen et al. 5 and Mordecki Q7|. Since P(M e(g) > k) = E[e-« T+ ( fc) ], the 
Laplace transform of the first passage time T + (fc), defined in (0 is equal to the 
left-hand side of l|33|) . 

The result below shows that it is possible to extend this result to a char- 
acterisation and description of the law of the up-crossing ladder process of X. 



Theorem 3 Let X be a Levy process of the form ^ such that —X^ ' is not a 
subordinator. Then the following are true: 

The ladder height process H is a subordinator given by 

N t 

H(t)=t + J2Un, (34) 

n=l 

where N is a Poisson process with intensity <5>(A*- + ^), where $(A ( - + ' ) ) is the largest 
root s of k X (-)(s) = A^ + ' and U n are i.i.d. random variables with distribution 

P(£4 e dy) = Vo ■ {S (dy), cMT {+) y)t {+) dy) (35) 
with So the delta measure in and rj given by fffij). 

Assume that the representation (m^ + \T^ + \ a^) is minimal. Then the cumu- 
lant k + o/(L _1 ,iJ) is given by 

K +( a , 3 )^ n ^ (+)(5 + ^ a)) (36) 
for s with 5ft(s) = 0, a > 0, where the denominator is taken to be one 
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Example 1 If there are no positive jumps (\( +s> = 0), and —X^ ' is not a 
subordinator, we find back (e.g. [Sj) that 

K + (a, s) = s + <&(a) 

where ^(a) is the unique positive real root of k X (-i(s) = a. 

From Theorems |3] and the extended form of the Wiener-Hopf factorisation 
(jl 1|> , we can now also determine the distribution of the downward ladder process 
by finding an explicit form for its cumulant. Inserting the explicit expression 
<|3(i|) for k + and comparing yields now: 

Corollary 1 Suppose the representation (m^, T^') is minimal. Then 
the cumulant kT of the dual ladder process (L, H ) is given by 

K {a,s)=C x(a-K{s) )= — . 

lliei<+)( s -PiW) 

where C > is some constant. 

Below we give elementary proof for Theorem 01 

Proof of Theorem We first determine the form of the process H. Since X^ 
has no positive jumps, AMt > implies that AXj: > and similarly, if 
dM t c > then dMj ' > as well, where denotes the running supremum 

of X^ up to time t. Write o\, 02, . . . for the jump times of X^ and define 
recursively for z = 1,2,... the stopping times Gi — inf{s > Di-\ : M s = X s } 
and Di = mf{erj : <jj > Gi}, where Go = 0, Do = cf\. Note that, for t € [Gi, Di), 
M t is continuous and that M may jump at Gi. As X^ is a compound Poisson 
process with rate A^ + ) , the differences Di — Gi are exponentially distributed with 
parameter A^ + ) and, as X^ has i.i.d. increments and M^'(e(M')) has an 
exponential distribution with mean $(A(+))^ 1 , M Di -M Gi are i.i.d. exponential 
with parameter $(A^ + '). Since we have taken the local time L to be equal to 
M c , the inter-arrival time of two jumps of H is exp(<i>(A < - + - ) )) distributed and 
thus the process H is given by 

Nt 

H(t) = M c {L-\t)) +X>M(L- 1 ( S ))l {AM(i - 1(s))>0} =t + J2 U n 

S<t 71=1 

where the U n are i.i.d. nonnegative random variables (since X is a Levy process) 
and N t an independent Poisson process with rate $(A^ + - ) ). The jump-size U n 
has the same distribution as the overshoot X T +( ) of X over the level 0, if 
X Q is distributed according to —A + B where A has the distribution £(dx) = 
P[(M(~) — -X^ _ ') e (A(+)) £ dac] 01 t ne distance of X^' to its supremum at an 
exponential time e(A^ + -*) and B, independent of A, is distributed according 
to the jump-sizes of X^ + \ Since the upward jumps of X are phase-type, it 
follows that the distribution of the overshoot X T ( Q \ has an atom in zero and on 
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(0, oo ) is (defective) phase-type (see e.g. Prop. 2] for a proof). A generator 
matrix of this phase- type distribution is seen to be given by T^ + \ with 'initial 
distribution over the phases' tj given by the distribution of the underlying 
Markov process at the moment of crossing. In the equivalent setting of the 
embedding (A, Y) = (Ax, Yx) of X, rj thus satisfies 

ijoCj) = P[r+ = j\A ~ £, Y ~ (0, q(+))] 



where j = 0, . . . , rrS + \ Thus the vector »7q is given by and l|35[) is proved. 

Finally, we turn to the proof of the identity (|36|) . On the one hand, since H 
is a compound Poisson process with unit drift, 

lim s _1 K + (a, s) = 1. 

s — >oo 

On the other hand, the form of the Wiener-Hopf factor y>+ and the fact that 
|X(+) | = m (+)+i imply that lim.^-oo s^+(s) = (-l)™ (+) n ieZ(+) ft/ det(T«). 
Combining with (|12|l completes the proof. □ 
We end this section with the characterisation of the form of the upcrossing 
ladder if X(~> is the negative of a subordinator. 

Theorem 4 Let X be a Levy process of the form 0) such that X^' is a sub- 
ordinator. Then the following are true: 

(a) H is a compound Poisson process with jump intensity c, the constant in the def- 
inition of L, and jump distribution of phase type with representation (t]q,T^ + ') 
where t]q is given by with (0,a^ + ') replaced by a( + \ 

(b) Suppose the representation [mS + \ a^,T^ + ') is minimal. For s with 5R(s) = 0, 
a > 0, 

.+ (,, -s) - c(l - ...l) 11 '^'" tff ' n det(r "' , v (37) 

det(sJ-T (+) ) ^^xwiPii )) 

where c > is the constant from the construction of the local time L. 

Proof of Theorem ^] Since in this case — X 1 -^ is a subordinator, X is irregular 
for (0, oo), that is T = T(0), the first entrance time of the positive half line, is 
positive a.s. and X will first enter (0, oo) by a jump. Thus, (T, Xt) are the first 
ladder time and ladder height respectively. The independence and homogeneity 
of the increments of X imply now that H increases by jumps which sizes are 
independent and distributed as Xt- By randomisation in the construction of 
local time L, we deduce that H is a compound Poisson process with intensity c 
and jump-sizes distributed as Xt- 

We now determine the law of Xt- Denote by Y the underlying Markov 
process of the jump AXt- Then Xt is the lifetime of Y with initial distribution 
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r] given in |J2BJ|. By the denning property of phase type distributions, Xt is 
has a distribution that is phase-type with representation (m^ , rj , T^). Since 
— X^' is a subordinator independent of X^ + \ it holds that a, the first jump 
time of X( + \ is exponentially distributed with mean 1/A^ + ) and 

/>oo 

V{X ( - ] e dx,Y a =k)= a fc A ( " ) e- A<_) *F(X 4 (_) 6 da;)di. 

Jo 

Conditioning on the position of at time <r and recalling that Y + , the 

upcrossing ladder process of Y, is a Markov process with generator Q + , it 
follows that 

poo 

Vo = « (+) j e«+"P(-X<-> +)) e dx), 

where as before e(q) denotes an independent exponential time with parameter 
q. Note that the equation for rj is the same as equation <|28[) (with a = 0) found 
before. 

Finally, we prove the identity i|37|) . Denote by G the distribution of T and 
note that in this case the Laplace transform G of G is E[e~ aT ] = r]' a l. By the 
randomisation in the construction of the local time L in this case the Laplace 
transform of L _1 (l) is seen to be 

oc 

K+(a,0) =E[c- aL ~ 1 ^] = P[Ni = n]G n — e^" 1 ), 

n=0 

where N denotes an independent Poisson process with rate c connected to this 
randomisation. Combining with the factorisation (|12[1 and the form l|33[l of the 
law of the supremum of X^ completes the proof in this case. □ 

5 Nonlinear iteration 

To solve explicitly for the first passage law, the law of the maximum l|33|) or 
the law l|35|) of the up-crossing ladder process of the Levy process X, we have 
to compute the generator matrix Q^- One approach achieving this proceeds 
by solving the equations (|2U[) and l|21|) numerically (e.g. via eigenvalue meth- 
ods). A different approach, exploiting the fact that in this case the matrix 
Qa has a special structure given by equations l|27[) - l|29|l . is to compute the 
sub-probability vector rj a . 

If is the negative of a subordinator, we consider the sequence (r)( n ') n 

where rj'°) a sub-probability vector and rf n+lS> , n > 0, is given by the right-hand 
side of (2EJ with the matrix Q + replaced by T (+) + i (+ V n) - If the distribution 
F(dx) — l^- a .<Qj.P(X ^ ^{x>o} aeTx tdx has exponential moments 

and non-zero mean, Corollaries 3.2 and 3.3 in Asmussen imply (rr')n con- 
verges geometrically fast to the solution of (|28|l . Asmussen [2] proves these 
results using a coupling argument. 
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In the sequel we therefore restrict to the case that X^' is not the negative 
of a subordinator. To prove the results we shall follow a route that is different 
from [2]. Set equal to the matrix in (|30|) and let S be equal to the set 

of sub-probability vectors in M m+1 with m = . From the interpretation as 
probability derived in the previous section it follows that the function tjj given 
by 

1> : V ^ ^J(T)( > a«)^<-)(-M« - m^ V ) (38) 

maps S to S. In the next result it is shown that the fixed point equation 
r] = ip(rj) is uniquely solved by the vector r/ Q . 

Theorem 5 Let a > 0. The following hold true: 

(i) The equation r) = ip(v) ^ as a unique solution r\ £ <S. 

(ii) With r?W = and rj^+V = ^(>7 (n) ), V (n) T Va asn^-oo. 

Proof (ii) Write r\ < r\' if tj' — rf is non-negative (coordinate- wise). We claim that 
ip satisfies the following monotonicity property: if r] < r\' then ip(i~i) < ip(r)'). 
Let P and P' be the transition matrices of the Markov chains with respective 
generators given by G(rf) :— M^ + ' + m^f] and G(rj'). Then the matrix 
D — P' — P satisfies the matrix differential equation 

D = DG(r]) + PW+> {n' - n) D(0) = O, 

the solution of which is given by D(t) = J Q P'(s)m < - + ' > (tj' — Tj)e^~ s ^ G ^ds. 
Hence, coordinate-wise, D is non-negative and since 

\(+) r 00 
m = ^p(0,a«) J e lG C?)p(-4; e dx) 

the claim follows. 

Starting with r/W = and setting t/™ +1 ) = ip(r]( n ') generates a sequence 
(v )n m <5 which is nonnegative, since the vector ?/>(0) is nonnegative, and 
coordinate-wise non-decreasing. Thus, the sequence has a limit in <S, say v, and 
by continuity of if) it follows that v = ip(v). 

(i) Since the matrix M^ + > + m^v, for v 6 5, is a generator matrix and 
solves equation i|19|) . uniqueness follows from Theorem [5] if k'(0 + ) > or a > 0. 
Indeed, in the latter case the generator matrix Q a is transient, whereas in the 
former case A satisfies sup t A t — +oo. 

If k'(Q + ) < 0, Theorem |21 does not apply and we need to provide a different 
argument to establish unicity. However, using exponential tilting, we shall show 
that this case can be reduced to the case of positive drift. If k'(0 + ) < 0, there 
exists a positive root 7 > of the equation k(s) = 0. Define the tilted measure 
P° for any positive constant c > 7 with k(c) < 00 by the Radon-Nikodym 
derivative 

dP c 

= cxp(cX t — K(c)t), t > 0, 



dP 
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and denote by k c , ip a ' c , <I> C and rf the respective quantities k, ip a , <& and r/ 
under the measure P c . It is straightforward to check that, for any c > 7, 
k c (s) = k{s + c), tp-' c {s) = tp-(s + c), $ c (a) = $(a) - c, k c '(0) = k'(c) > 0. 
Hence under P 7 the process X has a positive drift and unicity will follow if we 
show that l|31fl can be formulated in terms of quantities of the process X under 
the tilted measure P 7 . The next result (from [Tj) shows that, under the tilted 
measure P 7 , the jumps of X^ remain of phase type: 

Lemma 1 Under P 7 the jumps of X^ are of phase-type with representation 
given by 

(A (+ ' 7) ,a (+ ' 7) ,T (+ < 7) ) = (A (+) F (+) [-7],a (+) A/ J F , ( + )[- 7 ], A- 1 T (+) A + 7 J), 

where A is the diagonal matrix with kj on the diagonal where k = (7J — 
r (+))-i t (+). Further, = A'H^. 

Choosing the killing rate a 7 = a + A^ + '(l — 7]) and noting that a + A = 

a 7 + A 7 and recalling the form of from l|3U|) . we deduce that, under P 7 , 

the generator matrix of the upcrossing ladder process Y + is given by Q +r/ — 
M (+) ' 7 + m(+)< 7 r/ 7 , where 



M (+),7 = 



~$(a + A(+))+7 

A _1 tW A" 1 T (+) A + 7 / 



and mW' 7 = ($(a + A^) + 7,0) T , and r/ 7 satisfies 77 = V ,7 ( r 7) where ip 1 is 
defined as in (|38() but with a replaced by a 7 and all quantities by 'tilted' ones 
under the measure P 7 . Next we provide the link between ip and ip 1 . Writing 
Ai = diag(l, k\, . . . , & m (+>) and using Lemma^ it follows that 



\(+) 

(0, a (+ ))0- ( _, (-M (+) - mWj}) 



a + A(+) v ' J " x 



A (+)F(+)[- 7 ] rn „( +) ,^ ArU -, 7 ( r_ R/f (+)_ m (+) 
a + , 

A (+).7 



(( ; , (0, a (+) ' 7 ) ArV^l, (-MW - m(+)r7 - 7/) 



(0, a (+) ' 7 )</»-' 7 , (-A^(M(+' + 7 /)Ai - Ar^W^Ai) A7 1 



A (+)-7 



(0, a^ 7 )^:.' 7 , (-M (+) ' 7 - mW^JAj- 1 



a7 + A(+).7 v ' 

where r) = J/A^ Since k 7 '(0) > 0, Theorem |2 implies that rj — ijPijjj) has a 
unique solution fj in 5. Since it also holds that ip 7 (rj) is equal to i/j(t]) Ai and 
any solution of ip(r]) = r\ in S gives rise to a solution rj = rjAi of r) = ^(rj) in 
5, it follows that, also if «/(0) < 0, r\ a is the unique sub-probability vector that 
solves D 
If a > 0, it turns out that the convergence of the iteration to its solution is 
geometrically fast: 
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Theorem 6 Let a > 0. If n' ( _ ) (0 + ) > 0, the map ip is a contraction. If 



v x(-i 



(0 + ) < and k($(0)) < oo, the map ip^^°> is a contraction. 



Proof Writing \rj\ = XMI 7 ?*!}' A = \^ + \ a = and lo for the column vector 
that is one in state and zero else, it follows by the triangle inequality and the 
form of m(+) that for rj, rj' G S 



|V(r?)-V(V)l 
A 

< 



< 



< 



< 



a + A 
A 

a + A 
A 

a + A 
A 



(0,a) 



p(-4) eda; ) 



/•OO 

y (0, a)e^ M<+W+) «* W+) (»7 - 77') xPf-jfj G dx) 
y^ (0,a)e^ M(+) +" l<+) ^)l ia ;P(-/( (a ) ) e dx)*(a)\r, - r,'\ 



_ x n-l£ ) Ma)\r 1 -r,'\ 



for some vector £ x in the convex hull of r) and r;', where we used in the second 
line the mean value theorem and in the third line that the integrand is equal to 
x times the probability that some Markov chain is at time x in state 1. 



The proof is finished noting that E[-^ (a ) ) ]$(a) < 1 if and only if E[X[ 



> 



0. In case n' w-)(0 + ) < 0, the proof follows by replacing ip by ip®(°*> in above 



reasoning and recalling that k x <-)($(0)) = and /c*^)(0) > 



□ 
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